Magnetic dipolar ordering and quantum phase transition in Fe 8 molecular magnet. 
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We show that a crystal of mesoscopic Feg single molecule magnets is an experimental realization of 
the Quantum Ising Phase Transition (QIPT) model in a transverse field, with dipolar interactions. 
Quantum annealing has enabled us to explore the QIPT at thermodynamical equilibrium. The 
phase diagram and critical exponents we obtain are compared to expectations for the mean-field 
QIPT Universality class. 
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Quantum phase transitions 0,11 (QPT) have been ex- 
tensively studied in recent years. Physical realizations in- 
clude the superconductor insulator canonical phase tran- 
sition in cuprates [3|45|| , the onset of antiferromagnetism 
in heavy fermions [6l4lC|. t he p ressure driven insulator- 
metal transition in V2O3 [ll], and the order-disorder 



magnetic transitions driven by field (LiHoYF4 [l2j) or 
concentration (Cr K Vi_ K alloys [Hj]). In addition to their 
intrinsic interest, a plethora of new properties arise at 
nonzero temperature. 

In magnetism, the QIPT model is realized by a lattice 
of N coupled Ising spins in a transverse magnetic field 
H± , the Hamiltonian of which reads as follows 



N N 
i<j i 



(1) 



Here, a's are the Pauli spin operators, Jij the longitudi- 
nal couplings and As the ground-state tunnel splitting 
which depends on and vanishes with H±. The classi- 
cal long-range order that exists for H± = As — (say, 
ferro- or antiferro-magnetic) competes with the field- 
induced quantum fluctuations ([H,a] ^ 0). Long-range 
order is completely destroyed at the QIPT critical point 
when A S (H ± ) > A Sc where A Sc = 2JS 2 ~ k B T c and 
J = (l/Af) "Ar ^he ground-state becomes then a 
superposition of "up" and "down" spin states (here T c 
is the Curie temperature of a ferromagnet, which is the 
case for Fes). 

Crystals of interacting dipoles are natural candidates 
to observe a QPT 12]. Because of the weak interac- 
tion strength, the effects of quantum fluctuations can 



become noticeable at moderate transverse fields. In 
crystals of single- molecule magnets (SMMs), this effect 
is amplified because of very large inter-molecular dis- 
tances [T3 |. Furthermore, hyperfine interactions in 3d- 
based SMMs are generally weak enough not to block the 
QPT, whereas they play an important role in the case 
of the best known lanthanide-based insulators 15j . Fi- 
nally, SMMs show quantum phenomena like tunneling 
16l-|20[, interference [2 lj and superpositions of spin states 
22j . However, their complex mesoscopic character sug- 
gests that quantum spin fluctuations and therefore the 
QPT may, in these systems, be affected by decoherence 
[23| , therefore providing an interesting playground to in- 
vestigate QPT in the presence of a noisy environment 

A first indication of a possible QPT in a molecular 
crystal of Mni2 acetate SMM was inferred from mag- 
netic neutron diffraction experiments [25]. However, the 
experimental uncertainties close to M ~ hindered the 
exploration of the critical region. In addition, the local 
anisotropy axes of M1112 acetate molecules are slightly 
tilted due to the presence of several isomers giving rise, 
in the presence of a perpendicular field, to random fields. 
Very recently 26( , it has been argued that such fields turn 
Mni2 acetate into a realization of the classical random- 
field Ising model. The existence of a pure QPT in a 
crystal of molecular nanomagnets is, therefore, not yet 
established. Our aim here is to explore the experimental 
realization of the QIPT model in a disorder-free SMM 
crystal. 

We have chosen as a candidate the 
[(C 6 H 15 N3) 6 Fe 8 2 (OH) 12 ] 



as a 
system 



27] , hereafter 
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referred to as Feg. Each molecule has a net spin S = 10 
and a strong uniaxial magnetic anisotropy, which can be 
described by the following spin hamiltonian 



Hq — —DS Z + E (S^ — Sy 



(2) 



where D/k B = 0.294 K and E/k B = 0.046 K. Equation 
(|2"|) defines x, y and z as the hard, medium and easy 
magnetization axes. In the triclinic structure of Fe8, x, 
y and z axes are common to all molecules (28l |. 

Susceptibility experiments were performed on a single 
crystal [28j of approximate dimensions 3x2x1 mm 3 
and m = 4.72 mg. The complex magnetic susceptibility 
X = x'(T,u)) — ix"(T,u)) was measured down to 90 mK 
and in the frequency range 3 Hz < uj/2tt < 20 kHz using a 
homemade ac susceptometer thermally anchored to the 
mixing chamber of a 3 He- 4 He dilution refrigerator. A 
dc magnetic field H was applied with a 9 Txl Txl T 
superconducting vector magnet, which enables rotating 
i? with an accuracy better than 0.001°. 



The orientation of the magnetic anisotropy axes 28 1 
was first approximately determined through the identifi- 
cation of the crystallographic axes by X-ray diffraction. 
The magnetic easy axis z was oriented parallel to the 
ac excitation magnetic field of amplitude h ac = 0.01 Oe. 
The sample was completely covered by a non-magnetic 
epoxy to prevent it from moving under the action of the 
applied magnetic field. The alignment of H perpendicu- 
lar (±0.05°) to z and close (±20°) to the medium y axis 
was done at low temperatures (T = 2 K), using the strong 
dependence of %'(T, w) on the magnetic field orientation. 

Neutron diffraction experiments were performed on the 
thermal neutron two-axis diffractometer for single crys- 
tals D23 at the Institute Laue Langevin. A single crystal, 
thermally anchored to the mixing chamber of a a 3 He- 
4 He dilution refrigerator (T > 50 mK), was mounted so 
as to have the anisotropy axis orthogonal to the mag- 
netic field. Longitudinal field components arising from 
the nonperfect orientation of the crystal were compen- 
sated by means of a superconducting mini-coil inserted 
in the space between the dilution fridge and the super- 
conducting magnet. The intensities of different Bragg 
peaks, in particular (010) to which the magnetization 
components M z and M x strongly contribute with roughly 
equal weights and (—100) or (—212) to which M y and M x 
contribute significantly, were measured as a function of T 
(under a constant magnetic field) and H± (for a constant 
temperature) . The nonmagnetic contributions were mea- 
sured at high temperature (5 K) and H± = 0. The fit of 
the magnetic intensities enabled us to estimate the mag- 
netization components along the external magnetic field 
M± and along the anisotropy axis M z . 

The ac-susceptibility (Figure [IJ deviates from equilib- 
rium for low H± and low T, as shown by the vanishing 
of x' an d by the relatively large values of x" ■ Further- 
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FIG. 1: x (top) and x" (bottom) measured at u/2n = 333 
Hz and different H±. The top inset shows the shift of x" 
maxima with increasing H± . The solid line shows theoretical 
predictions for quantum spin-phonon relaxation that follow 
from Pauli's master equation as described in [3l|]. The bottom 
inset shows the magnetic relaxation time at /j.qH± = 1.7 T. 



more, the superparamagnetic blocking temperature Tb 
(associated with the x" maximum) is strongly frequency- 
dependent (not shown). In these dipolar Ising-like sys- 
tems, out of equilibrium conditions generally arise from 
the presence of energy barriers larger than the typical 
interaction energies (i.e. Tb > T c ), a situation which 
can be inverted by a quantum annealing process, via the 
application of a transverse magnetic field [2^, [5o|]. In- 
deed, when H± increases the height of the energy barrier 
decreases and spins are able to tunnel via lower lying 
energy levels, which leads to enhanced spin-lattice relax- 
ation rates. This explains qualitatively the decrease of Tb 
towards lower temperatures that we observe, as shown in 
the inset of Fig. [TJ 

At moderate H±, x! is n °t fully out of equilibrium 
(e.g. at fJ,oH± = 1.7 T where it remains finite down to 
lowest temperatures). The relaxation time r = x'Vx'^ es_ 
timated from these data is already very short and, more 
importantly, it remains constant below approximately 
250 mK. x' reaches full equilibrium for fioH± > 2 T 
as seen in Fig. [2] where x' measured at H± = 2.25 T 
becomes independent of frequency and temperature be- 
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low 340 mK (which we take as the critical temperature 
T c at this held). The saturation value corresponds to the 
inverse demagnetizing factor of our sample 1/N = 0.1 
emu/cm 3 Oe, as expected for an equilibrium ferromag- 
netic phase transition in which magnetic interfaces can 
move freely. Here, equilibrium arises because Ag is large 
enough (or/and the barrier small enough) to enable spins 
to relax via ground-state tunneling. 
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FIG. 2: x' measured at fioH± = 2.25 T. The data are in- 
dependent of frequency for any T > 0. The inset shows a 
log-log plot of demagnetization-corrected \\ vs the reduced 
temperature (for fj,oH± = 2.33 T with T c — 340 mK, o) and 
field (at T = 0.110 K with /i H c = 2.6 T, •). The slopes of 
the linear fits give the critical exponents 7 c i ~ 1.02(2) and 

7qu - 1-1(1) • 
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FIG. 3: H c -T c phase diagram. Solid symbols were obtained 
from the divergence of the equilibrium intrinsic susceptibility 
Xi- Open symbols were determined from the Curie- Weiss fit 
of 1/xi data measured in the paramagnetic region. The solid 
line represents a mean-field calculation of the phase boundary 
using Eq. @. The dashed line is the classical phase diagram, 
calculated from Monte Carlo simulations. Inset: x vs H± 
measured at T — 0.110 K and v = 333 Hz. 



observed Curie temperature of 0.57(5) K agrees very well 
with the Monte-Carlo prediction T c = 0.54 K [Hi ■ The 
solid line is obtained by solving numerically the mean- 
held Hamiltonian for the S — 10 spin 



We next show the existence of an equilibrium QPT 
and discuss its critical behavior. The intrinsic \[ cor ~ 
rected from demagnetizing effects is plotted vs reduced 
temperature, at 2.33 T, and held, at 110 mK, in the inset 
of Fig. O Under these conditions, \[ should follow, as it 
approximately does, the power laws 

where T c = 340 mK and [IqH c = 2.6 T are determined 
experimentally. The slopes give critical exponents 7 c i = 
1.02(2) and 7 qu = 1.1(1), in good agreement with 7 = 1 
of the Ising mean-held universality class in 3 dimensions. 
The H c -T c phase diagram shown in Fig [3] is obtained, 
above 2 T, by repeating the procedure described in Fig. 
H for different H±_ and, below 2 T, from the ht of (x0 _1 
to a Curie- Weiss law using the mean-held character of 
the phase transition (i.e. that 7d = 1 and that the Weiss 
temperature 9 ~ T c ). 

As expected, T c decreases when quantum fluctuations 
increase (with H±) showing a phase transition induced by 
thermal fluctuations at T C (H± = 0) = 0.57(5) K (semi- 
classical) and a phase transition induced by quantum 
fluctuations at /j,qH c (T = 0) ~ 2.6 T (quantum). The 



H = H - 9mHx(S x cos + S y sin0) - 2 J(S Z ) T S Z (4) 

where <f> is the angle between the magnetic held and 
the hard axis, 2 J is the molecular field coefficient and 
(S z )t is the thermal average of S z . We set J/fce = 
T C /2S 2 = 2.85 x 10" 3 K and = 70° that accounts 
for hqH c (T = 0) ~ 2.6 T. The agreement between the 
measured and calculated phase diagrams (Fig. [3]) is 
excellent despite that the held range where ferromag- 
netism is observed at equilibrium (when x' = can 
be very narrow as shown in the inset of Fig. |3l The 
dashed line represents the results of classical Monte Carlo 
simulations [331 ] which agree with classical mean-field 
predictions. The classical phase boundary is well ap- 
proximated by H C (T C ) = H c (0) [1 - T C /T C {H ± = 0)] 1/2 . 
In this model H c (0) equals the anisotropy held Hk = 
2[D-E (sin 2 - cos 2 0)] /gn B S ~ 3.8 T, which clearly 
overestimates the experimental critical field due to the 
absence of quantum fluctuations. 

Additional evidence supporting the existence of a 
phase transition to a dipolar ferromagnetic phase can be 
found in the results of neutron diffraction experiments. 
Figure |4] shows M 2 determined directly from magnetic 
diffraction intensities measured at noH± = 1.5 T. The 
inset shows held-dependent data measured at a fixed tem- 
perature of 0.14 K. Despite the high noise level, typical 
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of these experiments, the results show the onset of a net 
spontaneous M z in, respectively, the low-T (< 0.6(1) K) 
and low-i2j_ (< 3(1) T) regions. 
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FIG. 4: Squared longitudinal magnetization of a Fes crystal 
determined from neutron diffraction intensities measured un- 
der a constant fioH±_ — 1.5 T. The inset shows data measured 
at a fixed temperature T = 0.14 K as a function of H±. 

Summarizing, we have shown that the single molecular 
magnet Fes undergoes a dipolar ferromagnetic to param- 
agnetic phase-transition of (i) classical nature ((j,qH c = 
T and T c — 570 mK) and (ii) quantum nature (T c = 110 
mK and [ioH c = 2.6 T). Long-range magnetic order be- 
tween these mesoscopic SMM was enabled by the use of a 
quantum annealing procedure, which allows the spin sys- 
tem to attain thermal equilibrium despite the presence of 
high anisotropy energy barriers preventing spin-reversal. 
While the observed T c agrees very well with classical 
Monte Carlo calculations [32j], the H c — T c phase dia- 
gram shows important quantitative and qualitative de- 
viations from the expected classical behaviour at suffi- 
ciently low-T (or high-H±). By contrast, it is in excel- 
lent agreement with quantum predictions obtained via 
the diagonalisation of the mean-field Hamiltonian for the 
coupled S = 10 spins. The experiments show that the 
critical behavior and the phase diagram fulfill the expec- 
tations of the 3D Mean-Field Quantum Ising universality 
class. This agreement allows us to extend the notion of 
"macroscopic order parameter" 3J| of SMMs to the fer- 
romagnetic region. 
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